Non-gaussianity from the bispectrum in general multiple field inflation 
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We study the non-gaussianity from the bispectrum in multi-field inflation models with a general 
kinetic term. The models include the multi-field K-inflation and the multi-field Dirac-Born-Infeld 
(DBI) inflation as special cases. We find that, in general, the sound speeds for the adiabatic and 
entropy perturbations are different and they can be smaller than 1. Then the non-gaussianity can 
be enhanced. The multi-field DBI-inflation is shown to be a special case where both sound speeds 
are the same due to a special form of the kinetic term. We derive the exact second and third order 
actions including metric perturbations. In the small sound speed limit and at leading order in the 
slow-roll expansion, we derive the three point function for the curvature perturbation which depends 
on both adiabatic and entropy perturbations. The contribution from the entropy perturbations has 
a different momentum dependence if the sound speed for the entropy perturbations is different from 
the adiabatic one, which provides a possibility to distinguish the multi-field models from single 
£^ ■ field models. On the other hand, in the multi-field DBI case, the contribution from the entropy 

perturbations has the same momentum dependence as the pure adiabatic contributions and it only 
changes the amplitude of the three point function. This could help to ease the constraints on the 
DBI-inflation models. 

I. INTRODUCTION 

The inflationary scenario succeeds to explain the origin of temperature fluctuations of the Cosmic Microwave Back- 
ground (CMB). The increasing precision of the measurements of the CMB enables us to distinguish between many 
inflationary models. The primordial fluctuations generated during inflation are nearly scale invariant and gaussian. 
Thus the deviation from the exact scale invariance and gaussianity will give valuable information in discriminating 
many possible models. Especially, non-gaussianity of the primordial fluctuations will provide powerful ways to con- 
strain models (see e.g. [l[ for a review). The simplest single field inflation models predict that the non-gaussianity of 
the fluctuations should be very difficult to be detected even in the future experiments such as Planck [2| . If we detect 
large non-gaussianity, this means that the simplest model of inflation would be rejected. 

There are a few models where the primordial fluctuations generated during inflation have a large non-gaussianity. 
\Q • In the single field case, if the inflaton field has a non-trivial kinetic term, it is known that the non-gaussianity can be 
large. In if -inflation models where the kinetic term of the inflaton field is generic, the sound speed of the perturbations 
can be much smaller than 1 [1, HJ. This leads to a large non-gaussianity of the fluctuations. The Dirac-Born-Infeld 
(DBI) inflation, motivated by string theory, is another example [1, H, 0, IH ■ In DBI-inflation, the inflaton is identified 
with the position of a moving D3 brane whose dynamics is described by the DBI action. Again, due to the non-trivial 
form of the kinetic term, the sound speed can be smaller than 1 and the non-gaussianity becomes large [1, flo| . The 
third and fourth order actions for a single inflaton field with a generic kinetic term have been calculate d by properly 
taking into account metric perturbations and three and four point functions have been calculated fiol . [Til . Il2| . For 
the detailed observational consequences of single-field DBI-inflation see Refs. [ill EH, El, El, E3, El; El, HI] • 

Multi-field inflation models where the curvature perturbation is modified on large scales due to the entropy per- 
turbations 21] have been also extensively studied recently. In the case of the standard kinetic term, it is not easy to 
generate large non-gaussianity from multi- field dynamics [H, [23| (see however Rcf. [24]). In the DBI-inflation case, 
the position of the brane in each compact direction is described by a scalar field. Then DBI-inflation is naturally 
a multi- field model [2o| . The effect of the entropy perturbations in the inflationary models based on string theory 
constructions in a slightly different context is also considered in [H, [27| ■ Recently, Huang et. al. calculated the bis- 
pectrum of the perturbations in multi-field DBI-inflation with the assumption that the kinetic term depends only on 
X = —G IJ da4> I d^4) J l2 where </> J are the scalar fields (I = 1, 2, ...) and G/j is the metric in the field space as in the 
case of K-inflation [28]]. They found that in addition to the usual bispectrum of adiabatic perturbations, there exists 
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a new contribution coming from the entropy perturbations. Then they showed that the entropy field perturbations 
propagate with the speed of light and the contribution from the entropy perturbations is suppressed. This property 
can also be confirmed by the a naly sis of a more general class of multi-field models where the kinetic terms are given 
by arbitrary functions of X [2p j3Cj |. However, Langlois et.al. pointed out that their assumption cannot be justified for 
the multi- field DBI- inflation [3l|. Even though the action depends only on X in the homogeneous background, there 
exit other kind of terms which contribute only to inhomogeneous perturbations. They find that this dramatically 
changes the behavior of the entropy perturbations. In fact, it was shown that the entropy perturbations propagate 
with the same sound speed as the adiabatic perturbations. 

In this paper, we study a fairly general class of multi-field inflation models with a general kinetic term which includes 
K-inflation and DBI-inflation. We study the sound speeds of the adiabatic perturbations and entropy perturbations 
and clarify the difference between K-inflation and DBI-inflation. Then we calculate the third order action by properly 
taking into account the effect of gravity. Then three point functions at leading order in slow-roll and in the small 
sound speed limit are obtained. We can recover the results for K-inflation and DBI-inflation easily from this general 
result. 

The structure of the paper is as follows. In section II, we describe our model and derive the equations in the 
background. In section III, we study the perturbations using the ADM formalism. The second and third order 
actions are derived by properly taking into account the metric perturbations. Then we decompose the perturbations 
into adiabatic and entropy directions and write down the action in terms of the decomposed fields. In section IV, 
we study the sound speed in several models including K-inflation and DBI-inflation. It is shown that in general, 
adiabatic and entropy sound speeds are different and both can be smaller than 1 . It is shown that the DBI action can 
be obtained by requiring that the sound speeds for the adiabatic and entropy perturbations are the same. In section 
V, the third order action at leading order in slow-roll and in the small sound speed limit is obtained in terms of the 
decomposed fields. Then the three point functions are derived for a generalized model which includes K-inflation and 
DBI-inflation as particular cases. Section VI is devoted to the conclusion. 

II. THE MODEL 

We consider a very general class of models described by the following action 

^[M 2 Pl R + 2P(X IJ ,<f> 1 )] , (1) 
is the Planck mass that we will set to unity hereafter, R is the 

—g^d^d^, (2) 

is the kinetic term, t^,„ is the metric tensor. We label the fields' Lagrangian by P and we assume that it is a well 
behaved function. Greek indices run from to 3. Lower case Latin letters ■■■) denote spatial indices. Upper case 
Latin letters denote field indices. 

The Einstein field equations in this model are 

G>„ = Pg„ v + P^ud^d^ = T M „, (3) 

where Px 1 ' denotes the derivative of P with respect to X IJ [32j. The generalized Klein-Gordon equation reads 

g^(P x ud 1/ ct> I )^+Pj = 0, (4) 

where ; denotes covariant derivative with respect to <7 M „ and P j denotes the derivative of P with respect to (f> J . 

In the background, we are interested in flat, homogeneous and isotropic Friedmann-Robertson- Walker universes 
described by the line element 

ds 2 = -dt 2 + a 2 {t)5 l] dx 1 dx j , (5) 
where a(t) is the scale factor. The Friedmann equation and the continuity equation read 

3H 2 = E , (6) 



d A xJ~- 



where cf) 1 are the scalar fields (/ = 1,2, ...,N), Mpi 
Ricci scalar and 



X IJ = 



3 



E = -3H (E + P ) , (7) 
where the Hubble rate is H = a/a, Eq is the total energy of the fields and it is given by 

E a = 2X* J P XIJ - P , (8) 

where a subscript zero denotes background quantities and Xq J = 1/20q0q . The equations of motion for the scalar 
fields reduce to 

Pxuft + (zHPxu + V) - P. j = 0. (9) 
III. PERTURBATIONS 

In this section, we will consider perturbations of the background ([5]) beyond linear order. For this purpose, we 
will construct the action at second and third order in the perturbations and it is convenient to use the ADM metric 
formalism [H M, M, H, H, [H . The ADM line element reads 

ds 2 = -N 2 dt 2 + hij (dx l + N l dt) (dx 3 + N j dt) , (10) 

where N is the lapse function, N l is the shift vector and hij is the 3D metric. The action (jTJ) becomes 

S=\ J dtd 3 xVhN ( (3) i? + 2P(X IJ , (j) 1 )^ + i J dtd^VhN- 1 (EijE ij - E 2 ) . (11) 



The tensor Eij is defined as 



Eij = i (hij - ViNj - S7jNi ) , (12) 



and it is related to the extrinsic curvature by Kij = N 1 Eij. V, is the covariant derivative with respect to hij. X IJ 
can be written as 

1 AT -2 

X IJ = -h";),o'i),o J + —v'v 1 , (13) 

where v 1 is defined as 

v 1 EE ft - NiVjcf) 1 . (14) 
The Hamiltonian and momentum constraints are respectively 

(3) i? + 2P- 2N- 2 P iX uv I v J - N- 2 (E VJ E lj - E 2 ) = 0, 

Vj (n^EI^ - V, (N- l E) = N-^Pxvv^i^. (15) 

We decompose the shift vector N l into scalar and intrinsic vector parts as 

Ni = N t + ditp, (16) 
where diN -1 = 0, here and in the rest of the section indices are raised with 5y . 

A. Perturbations in the uniform curvature gauge 



In the uniform curvature gauge, the 3D metric takes the form 

hij a Sij , 

<j> I {x,t) = ^{t)+Q I {x,t), (17) 



where Q 1 denotes the field perturbations. In the following, we will usually drop the subscript "0" on 0q and simply 
identify tfi 1 as the homogeneous background fields unless otherwise stated. 
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We expand N and N l in powers of the perturbation Q 1 

N = l + ai +a 2 + -- 
Ni = Ni™ + N z {2) + 



if) = ipt + ip 2 + 



(18) 

(19) 
(20) 



~ (n) 

where a n , Ni and ip n are of order (Q 1 ) 71 . At first order in Q 7 , a particular solution for equations ([15)) is: 



5 2 Vi = 



2# 



6H 2 ai 



+ P,kQ K - <P I ^ J P x u K Q K + (P x w + ^^P x «x«) - ^Q 1 )] • (21) 



The second order action is calculated as 

,3 



'(2) 



, or 
dtd x— 



3 : , 



X[ j X[ m P x u x lm + P.xuQ'Q' 1 - a- 2 P, XI .,d l Q I d t Q J - -^^Px'-'P.x^'Q'Q™ 



H 



'P xL uXf M + P, K Q K ) + 2P tXI ., K Q K (-^0 J ai + tfQA + P : klQ K Q 1 



+P }X u Utf^al ~ 4ai^Q J ) 



(22) 



where 



X[ J = -a^ 1 ^ + (/ Q J) . 



(23) 

After integrating by parts in the action and employing the background field equations, the second order action can 
be finally written in the rather simple form 

1 



5(2) = 



dtd 3 xa 3 [(P, x u + P^iKxjL^^QtQ 1 



1 



,P. x ,.,d l Q I d l Q J - MuQ T Q J +N I jQ J Q I ] , 
with the effective squared mass matrix 

M u = -P.U + 



(24) 



X LM 



H 



{P,x JK P.ix LM + P,x ik P.jx lm ) 



--^x MN x p ®p x m NxPQ p x j K p x j L 4> k j> L 



1 d 
~a?~db 



H 



-P x IK P x JL( 



'kKIL 



-MN 



X 

NlJ = 1\Pj X IK — — P X IK X MN P,X JL< > 



H 



(25) 
(26) 



In the same way, the third order action is given by 

'(3) 



Stx\ = I dtd 3 xa 3 



3H 2 a 2 + ^ai^Vi + ^ (dVi^Vi - did^d^x) 



Oil 



- \a 3 tf$ J + aU'Q J + a-^axftd^xdiQ 3 - ^-Q J Q J - a~ 2 d t Q J ( Q z 0tyi + ^d'Q 1 



1 



.2 11 1 J 



,.2111 J 



1 



1 



1 ArA 



1 



vLM r> 

A l r_ x iJX Ll ' 

P.x'->kQ K 



-—PljLJ y + -A x A x Aj 1 j x ij x lm x qr + -A 1 A x ^ x i'X LM KW 

lb I 
~X[ J P x u LM Q L Q M + \P, IJK Q I Q J Q K 



(27) 
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B. Decomposition into adiabatic and entropy perturbations 



We can decompose the perturbations into the instantaneous adiabatic and entropy perturbations, where the adi- 
abatic direction corresponds to the direction of the background fields' evolution while the entropy directions are 
orthogonal to this [37]. For this purpose, following [2{|, we introduce an orthogonal basis e^(n = 1,2, ...,N) in the 
field space. The orthonormal condition is defined as 

P x ue I n e m = S nm , (28) 
so that the gradient term P x udiQ I d' l Q J is diagonalized [38|. We pick up the adiabatic vector as 

e[ = ^ , ; , (29) 

1 P X .,K(j) J (t> K 



which satisfies the normalization given by Eq. (|28p . The field perturbations are decomposed on this basis as 

Q 1 = Q n ei . (30) 
Wc defined the matrix Z mn which describes the time variation of the basis as 

= e m^rnn , (31) 

which satisfies Z mn = —Z nm — P x ue m e J n as a consequence of (Px /je n e m) = 0- 
In terms of the decomposed fields, the second order action (|24|) can be rewritten as 



If 1 

S(2) = 2 dtd 3 xa 3 [K, mn (D t Q m )(D t Q n ) ^8 mn diQ m d l Q n 

-M in 1 1 QmQn + NmnQn{DtQm)\ , (32) 

where 

DtQrn — Qm ~\~ Z mn Q n , (33) 

K mn = 8 mn + {P^ X M N ^ M 4> N )P^ XI K X ., L e[e^e(e^ 1 (34) 
M mn = M u e m e J n , (35) 
M mn = M u e m e J n . (36) 

From the constructions, JC mn , -Mmn and Mmn are symmetric with respect to m and n. The explicit form of the 
effective squared mass matrix in this representation is 

M mn = + 4 {P. X KL4> K ^ L f/ 2 {P, mX M N e^e^)5 nl 

rl 

1 r (P x kl j> K j> L f{P tX MM X P Q ef ef ef ef)6 ml 6 nl 



AH 2 
1 d 



e m e n ! (37) 

III in 

1 



- (P x mn <P M cf>" )P x <k P XJL ef e L x 



M mn = --(P x p Q ^ P (l> Q ) 2 {P X KL X MNe^efe^)S nl 
Id 



+2JP tXLM <}> L 4> M {P, n x^e m ef), (38) 



where P mn = Pjje l m e J n and P nX iK = P Jx i K e J n . 
The equation of motion is obtained as 



^g-— [a 3 (2K. mr DtQ m + A/mrQm)] — (2/Cmn-^nr + N mr )Dt( 



(2M mr + MmnZmr) Qm + Qv = . (39) 

or 
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IV. LINEAR PERTURBATIONS 

In this section, we study the linear order perturbations using the second order action derived in the previous section. 

A. K-inflation 

Let us consider K-inflation models where P(X IJ , (j) 1 ) is a function of only the trace X = X IJ Gu(4> K ) of the kinetic 
terms where Gi,j{<p K ) is a metric in the field space: 

P{X IJ ,tf)=P{X,cj> 1 ). (40) 

The derivatives of P can be evaluated as 

P x u = GuP x , (41) 

Pi = Igjkj^^Px+Pj, (42) 

P,x lJ x KL — GijGklP,xx , (43) 

P X "k = \glm,k^ L ^ M GijP,xx + G LlK P x + G u Pxk , (44) 



1 



Pjj = ^Gkl.iGmnj^^^^Pxx + ^Gkl,iA K ^ L P,x 



+ -<t> M iV (Gmn,jPxi + GmnjPxj) + Pu , (45) 



and the sound speed is defined as 



2 Px 



P x + 2XP XX 

In terms of the decomposed field, the second order action can be written as 

1 



(46) 



s (2} = ; / dtd 3 - 3 



\ - 1 ) 8 lm Sln \ (D t Qm)(D t Q n ) 



(47) 



where we do not show the explicit forms of M mn and Af m n- 

The sound speed agrees with the adiabatic sound speed defined by c 2 s = dP/dE. The fact that the sound speeds for 
the entropy perturbations are unity has been recognized in Ref. [29] . This is because the non-trivial second derivative 
of P only affects the adiabatic perturbations, which is the consequence of the fact that the entropy field s satisfies 
s = in the background and thus it has no first order perturbations. 

B. DBI-inflation 

An interesting class of models is the DBI-inflation which describes the motion of a D3 brane in a higher dimensional 
spacetime. The DBI action is given by 



S = - 



d^j^-yj-detig^ + fi^Gui^d^du^]. (48) 



Recently it was pointed out by Ref. [3l| that the multi-field DBI-inflation is not included in the multi-field K-inflation 
discussed in the previous subsection. Indeed, P(X IJ ) is not a function of X, but it is given by 

P{X IJ ,<j ) I )=P{X,<p I ) : X=^p-, (49) 
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where 

V = det{8 I j-2fX I J ) 

= 1 - 2fG u X IJ + 4f 2 X l /Xj ] - 8fX l /xjXx ] + m,f 4 X l /xjX^X^ ] . (50) 

In the background, X = X. However, this does not mean that the full action is a function of X only. The DBI action 
takes a specific form of P 



f 

where we allow for a potential V^ 1 ). The sound speed is defined as 



P(x, 4, 1 ) = -i ( Jl - 2fX - 1 ) - V& 1 ), (51) 



c 2 s = -= . (52) 

P.x + 2XP ,XX 



The derivatives of P can be calculated as 



- pjt y-^n ) . (53) 



P *" XKL = P *4^)(j^) +P 4 dX^fxKL > (54) 



where 



dX 



dxu - clGu + 2/Xu, (55) 
2/ (gjjGkl - \g ik G jl - l -G IL G JK ] + 0(X /J ) . (56) 



dX IJ dX KL " 2 14 " 2 

Here we do not explicitly write down the higher order terms in X 7,7 in the second derivative as they will not contribute 
to the final result. In the following, we will omit these terms. We can also show that 



Pi = -G JKJ <j> J <j> K P,x+Pi> (57) 
Pjj = \g K ljGmn^ K ^ L ^ M ^ N 'P xx + \g KL jj^ k 4> l P x 

1 



+ -</> A V v {G M n,jP iXI + GmnjP^j) + Pjj , (58) 



\^MN,jr.xi "•" ^ MN < IJr ,Xj) ^VJ 

P.x ij k — ((1 — 2/X)Gjj ; js: — 2/^kXGu — 2JGlm,kX lm Gjj 
1f,KXu + 2JGjl^kX L j + 2fGjM,KXj M )P x 



+(c 2 s Gu + 2fXu) ( -G LM ,Ki L i M P xx + P.xk ) • (59) 



1 

-( 

2 

It is worth noting that even though P t x IJ K seems to be a bit complicated, we can show that 

P,x"k¥ = \glm,k4> L ^ M 4>iP^xx + frPxK + Gij, K j> J P,x > (60) 

which is just the same form as the K-inflation case. We can also show that 

F x ,.,rt J = 2JP 1 . (61) 

The orthonormality conditions for the basis give 

1 1 1 - c 2 
eie m i = -~ -S mn - -g j^mi^i • ( 62 ) 
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Using these results, the second order action can be written in terms of the decomposed perturbations as 

'(2) 



.S'.o, = - / <ll (I'm'' 



—oS mn (D t Q,n)(D t Q n ) T:8 mn diQ m d l Q r . 



(63) 



Unlike the K- inflation models, all field perturbations have the same sound speeds as was pointed out by Ref. [3 11 ] - 
In order to understand the difference between K-inflation and the DBI-inflation, we will consider a generalized model 
where both cases are included. 

C. Generalized case 

Let us consider models described by 

P{X IJ ,cj> I )=P{Y,<p I ), (64) 

where 

Y = G u (<t> K )X IJ + ^(X 2 - Xfxj) . (65) 

The functional form of Y is chosen so that Y = X = GijX ij in the background as in the DBI-inflation model. This 
model includes as particular cases the K-inflation model for b = and the DBI-inflation for b = —2/ and if P has 
the DBI form. This might be surprising as the DBI action contains additional terms of order f 2 and / 3 in X (see 
equations (|49|) and ([50)) '). but it turns out that these terms do not contribute to the second order action and the 
leading order third order action. 

Following a similar procedure to the previous subsection, the second order action can be written in terms of the 
decomposed perturbations as 



-5(2) = - / (Ihf'.r/r' 



) > (D t Q m )(D t Q n ) 

ry 1 + OA 



+ NmnQn(DtQm) 



a 2 



(66) 



Now we are in a position to explain the difference between K-inflation and DBI-inflation. As in the K-inflation case, 
the non-trivial second derivative of P affects only the adiabatic perturbations. On the other hand, the non-linear 
terms of X IJ in Y only affects the entropy perturbations as they vanish in the background. Then the sound speed 
for adiabatic perturbations c 2 ad and for entropy perturbations (? en are given by 

c 2 ad EE -= , <4 = 1 + bX , (67) 

Fy + 2XP <YY 

and they are independently determined by P^yy and d 2 Y / (dX 1 J dX KL ) respectively. Thus in general they are different. 
Let us derive the condition under which the two sound speeds are the same, i.e., c 2 d = c 2 n . This condition is given by 

Pyy bX , . 

2X^- = — . 68) 

P y 1 + bX y ' 

Then we find that the DBI action is a solution of this equation where b = —2/ [39]. 



V. THE LEADING ORDER IN SLOW-ROLL THREE POINT FUNCTION 



In this section, we will calculate the leading order in slow-roll third order action for the generalized model of the 
previous subsection and then we shall calculate the leading order three point function for both adiabatic and entropy 
directions. Finally we will obtain the three point function of the comoving curvature perturbation. 
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A. Approximations: slow-roll 



In order to control the calculations and to obtain analytical results we need to make use of some approximations. 
We will use the slow-roll approximation, where we define a set of parameters and assume that these parameters are 
always small until the end of inflation. We define the slow-roll parameters as 

- & XPy - k (f\Q\ 

Cad C en /nn\ 

Xad= 77, Xen= 77 ■ (70) 

C-adtl C erl -tl 

It is important to note that these slow-roll parameters are more general than the usual slow-roll parameters and that 
their smallness does not necessarily imply that the fields are rolling slowly. Assuming that the parameters Xad an d 
Xen are small implies that the rates of change of the adiabatic and entropy sound speeds are small, but the sound 
speeds themselves can have any value between zero and one. 

It is convenient to define a parameter that describes the non-linear dependence of the lagrangian on the kinetic 
term as 

X=^X 3 P tYYY + X 2 P. YY . (71) 
We will also assume that the rate of change of this new parameter is small, as given by 

< 72 > 

At the end of this section, we will show that the size of the leading order three point function of the fields is fully 
determined by five parameters evaluated at horizon crossing: e, A, H and both sound speeds. 

It turns out that the equations of motion for both adiabatic and entropy perturbations at first order form a coupled 
system of second order linear differential equations, see Appendix for details. In general, the coupling (denoted by 
£ in equation (|A10[) ) between adiabatic and entropy modes cannot be neglected but in this work we will study the 
simpler decoupled case, where we assume that £ is small when the scales of interest cross outside the sound horizons, 
i.e., we will assume that £ ~ 0(e)- With these approximations the adiabatic and entropy modes are decoupled and 
the system of equations of motion can be solved analytically. For simplicity, we will also assume that the mass term 
present in the entropy equation of motion is small, i.e., [i 2 /H 2 <C 1 (refer to Appendix for more details). When 
calculating the leading order three point functions, we assume that the quantities related to the time derivatives of 
the basis vectors given by Z mn are also slow-roll suppressed. Finally, the calculation of the three point functions in 
the next subsections is valid in the limit of small sound speeds. Our results will also include sub- leading terms of 
0(1) but these terms will in general (for small sound speeds) receive corrections coming from terms of the order of 
e/Cg, that we have neglected. 

B. Third order action at leading order 

At leading order in the previous approximations, the third order action for the general model (JTJ is calculated as 

o 

\p x ik x jl^ j r Q K ^diQ- J d l Q L . (73) 



#(3) = \ I dx 3 dta 3 



2 



After decomposition into the new adiabatic/entropy basis the third order action can be written as 



5(3) = / dx 3 dta 3 



2 2a z 



(74) 
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where we define the coefficients ^ nm i and T„ m ; as 

2„ m , = P XI K X j L ^P^4^e[ I efJeief ) 

+lp^*x™x*»iPx^P)*<M X ep#e%4»e?>, (75) 

T„ m , = P, XI K X j L y/p X M N 4> M j> N e[ I eVeie{' . (76) 

We shall now give some useful formulae of the previous quantities for the different inflationary models considered in 
this work. 

1. K-mflation 

For the K-inflation model we have 

" fnvrt \ — -( 2XP X X <- j- . 4: X 2 P XXX c £• r 1 , n n\ 

^nmi = (2XP.x) 2 — ~ o 1 ^ n d m i) + - ~ o nl d ml dn , (77) 

V "< x ) 

T mnl = (2XP x )-i 2XP ' XX S nl 6 ml . (78) 
P,x 



2. DBI-inflation 



For the DBI-inflation scenario they are given by 

il-C? 



E nml = (2XP x )-^—^S 1(n 5 ml) , (79) 

_ L / 1 — c 2 1 - C 2 \ 

Tnmi = (2XP X ) 2 I ^-^nl^mi - 2 ^ (£„l<iW - 5 n (m^)l) ) ' ( 80 ) 

where c 2 s should be understood as the sound speed defined in Eq. (1521) . 



3. Generalized case 

For the generalized case of subsection II V Cl S nm ; and T nm ; can be written as 



(2IP y )-5 



o 2 di ( „d TOi) + ~ -3—5 dnldmldU 

C ad C en V *,y C ad C en / 



„2 _ J _ J2 



(81) 



(2XP y )-i ( l^«5 nl( 5 m; - 2(1 2 C ' J - S n{m 8 in ) ) , (82) 



and it is obvious that the DBI-inflation is a specific case of the general model with c ad = c en = cj 

C. The three point functions of the fields 

In this subsection, we derive the three point functions of the adiabatic and entropy fields in the generalized case and 
at leading order in slow-roll and in the small sound speeds limit. We consider the two-field case with the adiabatic 
field a and the entropy field s. 

The perturbations are promoted to quantum operators like 

Q»(r,x) = J d 3 kQ„(r,k)e* kx , (83) 
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where 

Q„(r,k) = u„(r,k)a n (k) +<(r,-k)at(_ k ). (84) 
a„(k) and (— k) are the annihilation and creation operator respectively, that satisfy the usual commutation relations 
[MkO^LM = (2vr) 3 ^ 3 )(k 1 -k 2 )S nm , [a n (k 1 ),a m (k 2 )] = [4 ( kl ), <4(k a )] = 0. (85) 
At leading order the solution for the mode functions is (see Appendix for details) 

u„(r,k) = A n ^(l + tkc n T)e' lkc " T , (86) 

where c n stands for either the adiabatic or the entropy sound speeds. 
The two point correlation function is 

(0|Q„(r = 0,kx)Q m (T = 0,k 2 )|0) = (2^) 3 ,5 (3) (k 1 + k 2 )V Qn ^-S nm , (87) 

where the power spectrum is defined as 

and it should be evaluated at the time of horizon crossing c„»fci = a*H* 

The vacuum expectation value of the three point operator in the interaction picture (at first order) is [34|, |4^| 

(n|Q,(<,k 1 )Q TO (t,k a )Q n (t,ka)|fi) = -i / dt(0\ [Q,(t 1 ki)Q TO (t,k a )Q n (t,k 5 ),.Bj(t)] |0>, (89) 

where to is some early time during inflation when the field's vacuum fluctuation are deep inside the horizons, t is 
some time after horizon exit. |0) is the interacting vacuum which is different from the free theory vacuum |0). If 
one uses conformal time, it's a good approximation to perform the integration from — oo to because r rs — (a£f) _1 . 
Hi denotes the interaction hamiltonian and it is given by Hi = — L^, where L 3 is the lagrangian obtained from the 
action (174]) . 

At this order, the only non-zero three point functions are 

(Q|Q ff (o,k 1 )g CT (o,k 2 )g CT (o,k 3 )|n) - (2^) 3 <5( 3 >( kl + k 2 + k 3 ) 2c °^' 6 ^ 1 

M k l k i n 7.2,, f-, , k 2 + h . n k 2 kt 



6&(C + G 4 )^ - C lfc ?k 2 • k 3 ( 1 + *p + 2^ 
+2 cyclic terms 



(90) 



(fi|Q CT (0,^)0,(0, k 2 )Q s (0,k 3 )|O) = (27r) 3 ( 5( 3 )(k 1 +k 2 + k 3 )^4^ 



nf =1 fc 3 £ 

en a 1 "2 

C 2 c en fc|ki • k 2 ( 1 



K K 2 

„2 ,2i 



2(d+C 2 )«k 3 .k 3 1 



+ (&2 fa) 

+ 2c^ 



k 2 + h , „ 2 fafa 



A" 



A 2 



(91) 



where A = fci + fc 2 + fa, A = c ac ;fci + c e „(fc 2 + fa), cyclic terms means cyclic permutations of the three wave vectors 
and (fc 2 fa) denotes a term like the preceding one but with fc 2 and k 3 interchanged. The pure adiabatic three point 
function is evaluated at the moment r* at which the total wave number A exits the horizon, i.e., when Ac 0( j* = a*-H*. 
Because of the different propagation speeds, the adiabatic and entropy modes become classical at different times, 
however at leading order we assume that the background dependent coefficients of (|9ip do not vary with time and so 
they can also be evaluated at the moment r*. 



12 



The different constants Cn are given by 



d = {2H 2 e)~* 1 f ad , C 2 = -2{2H 2 e)-? 1 ^ 



c 



2 ' ^ -v,"" v., 2 

ad °en 



C ad C en \ e C ad C en 



D. The three point function of the comoving curvature perturbation 

In this subsection, we calculate the leading order in slow-roll three point function of the comoving curvature 
perturbation in terms of three point function of the fields obtained in the previous subsection. 
During the inflationary era the comoving curvature perturbation is given by 

n= HQ^ 



It is convenient to define the entropy perturbation S as 




so that Vs, — V-r, , where the subscript * means that the quantity should be evaluated at horizon crossing. 

In this work we will ignore the possibility that the entropy perturbations during inflation can lead to primordial 
entropy perturbations that could be observable in the CMB. But we shall consider the effect of entropy perturbations 
on the final curvature perturbation. We will follow the analysis of Wands et. al. [4l[ , where it has been shown that 
even on large scales the curvature perturbation can change in time because of the presence of entropy perturbations. 
The way the entropy perturbations are converted to curvature perturbations is model dependent but it was shown 
that this model dependence can be parameterized by a transfer coefficient Tns SH like 

K=R,+%5 t =AQ„+i s Q s », (95) 

with 

A, = I — %= ) . A, - - T RS I — =,/— I • (%) 




Using the previous expressions we can now relate the three point function of the curvature perturbation to the 
three point functions of the fields obtained in the previous subsection. The three point function of the curvature 
perturbation is given by 

(^(kOftCk^ka)) =^(Q £r (k 1 )Q CT (k 2 )Q CT (k3))+^ CT ^({Q £r (k 1 )Q s (k 2 )Q s (k 3 ))+2perms.). (97) 

For the DBI-inflation case the previous equation can be simplified and the total momentum dependence of the three 
point function of the comoving curvature perturbation is the same as in single field DBI 31]. For our general model 
this is no longer the case, i.e., the different terms of the previous equation have different momentum dependence. Once 
again one can see that DBI-inflation is a very particular case and more importantly it provides a distinct signature 
that enables us to distinguish it from other more general models. 



VI. CONCLUSION 



In this paper, we studied the non-gaussianity from the bispectrum in general multi-field inflation models with a 
generic kinetic term. Our model is fairly generic including the K-inflation and the DBI-inflation as special cases. We 
derived the second and third order actions for the perturbations including the effect of gravity. The second order 
action is written in terms of adiabatic and entropy perturbations. It was shown that the sound speeds for these 
perturbations are in general different. In the K-inflation the entropy perturbations propagate at the speed of light. 



13 



The DBI-inflation is a special case where the sound speed for the entropy perturbations is the same as the adiabatic 
sound speed. We found that, from the requirement that the sound speeds for adiabatic and entropy perturbations are 
the same, we obtain the DBI form for the action. 

Then we derive the three point function in the small sound speeds limit at leading order in slow-roll expansion. In 
these approximations there exists a three point function between adiabatic perturbations Q a and entropy perturbations 
Qs, (Q a {k.i) Q s O^z) Qs {^3)), m addition to the pure adiabatic three point function. This mixed contribution has a 
different momentum dependence if the sound speed for the entropy perturbations is different from the adiabatic one. 
This provides a possibility to distinguish between the multi-field models and the single field models. Unfortunately, 
in the multi-field DBI case, the sound speed for the entropy perturbation is the same as the adiabatic one and the 
mixed contribution only changes the amplitude of the three point function. This could help to ease the constraints 
on DBI-inflation as is discussed in Ref. [3l| . 

In order to calculate the effect of the entropy perturbations on the curvature perturbation at the recombination, we 
need to specify a model that describes how the entropy perturbations are converted to the curvature perturbations. 
In addition, even during inflation, if the trajectory in field space changes non-trivially, the entropy perturbations can 
be converted to the curvature perturbation. In this paper, we assumed that this does not happen and neglected a 
mixing between the curvature and entropy perturbations. It would be interesting to study this mixing in specific 
string theory motivated models. 

Note added: While we were writing up this work, similar results appeared on the arXiv 44 1 . 
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APPENDIX A: EQUATIONS OF MOTION FOR THE FLUCTUATIONS 



Here, we derive the equations of motion for linear perturbations for the generalized model introduced in IIV CI In 
terms of the field space "covariant quantities" |42j which are given by 



II — 21 



U IK 



VtQ 1 ee Q 1 + Y 1 JK 4> J Q K 



VjVjP EE Pjj - TfjP iK , 

■T~) I T^Z T^I I T^I T^Af T^I A'l 

/c KLJ = 1 KJ,L ~ 1 KL,J ~r 1 LM 1 JK ~ 1 JM L LK ; 



(Al) 
(A2) 
(A3) 
(A4) 



(Tj K denotes the Christoffel symbols associated with the field space metric Gjj), the second order action can be 
expressed as 

1 



5(2) = - I dtd i xa i 



J dtd 3 xa 3 (PyG/j + P,YYhh) V t Q I V t Q-' 



1 ^ 



—P,y [(1 + bX)G u - bXu] d t Q J d l Q J - M U Q J Q J + 2P, YJ( j )I Q J V t Q 1 



(A5) 



with the effective squared mass matrix 
M u -- 



VjVjP - p, Y n IK Lj4> K ^ L 

Is-,; 1 



XPv ~ 



XP 3 Y 

' 2H 2 



(1 











c ad/ 


<I>i4>j 



(A6) 



It is worth noting that except for the coefficients of the kinetic term and the gradient term, this action is the same 
as the K- inflation case and DBI-inflation case which are derived in [29( and [3l|, respectively. 



14 



From now on we will derive the equations of motion for the fluctuations. For simplicity, let us now restrict our 
attention to the two field case (I = 1, 2). Then, the perturbations can be decomposed into Q 1 — Q^e^ + Q s e a , where 

= e{ and el is the unit vector orthogonal to e£. As in standard inflation, it is more convenient to use conformal 
time t = J dt/a(t) and define the canonically normalized fields [43| 

a 



—Qo 

Cad 



— Q s . 

Cert 



(A7) 



From the similar calculations with K-inflation and DBl-inflation cases analyzed by [29| and [311 ] . respectively, we 
find the equations of motion for v a and v s as 



,2 k 2 
'ad K 



Z 
Z 



v. 



Cen k \-a (J, S )V S £v a 

a 



z 
z 



0. 
0, 



where the primes denote the derivative with respect to r and 

a 



oP,YCad 
p 

a a 



(1 + c 2 ad )P s - c 2 ad a 2 P ys 



-g 2 R 



CadH 



1 P's 
Kd* P 2 y 

t = a\ I Py ! 



P.YsP.s 
' P 2 



with 



(A8) 
(A9) 

(A10) 
(All) 
(A12) 



a = V2X , 



Ps = Pjei 



'P.YC 



p 



Ys 



P 



,YI e s 



'P.YC, 



P ss = (V I -DjP)e I s e J s P, Y 



(A13) 



and R denotes the Riemann scalar curvature of the field space. 

If we assume that the effect of the coupling £ can be neglected when the scales of interest cross the sound horizons 
the two degrees of freedom are decoupled and the system can be easily quantized. If we further assume the slow-roll 
approximations, the time evolution of H, c a d, and a is small with respect to that of the scale factor and the relations 
z"/z ~ 2/r 2 and a" /a ~ 2/t 2 hold (see section V.A for these approximations). The solutions of (|A8[) and (|A9[) with 
the Bunch-Davies vacuum initial conditions are thus given by 



1 



Vak 



Vsk 



V^kCad 
1 



-ikc a dT 



— ike 



when /j, 2 / H 2 is negligible for the entropy mode. 

Therefore, the power spectra for Q a and Q s are obtained as 



1 - 
1 - 



kc ad T 
i 



(A14) 
(A15) 



H 2 



^ 2 C a d 



H 2 



An 2 c e 



(A16) 



which are evaluated at sound horizon crossing. The ratio of the power spectra for the adiabatic and entropy modes 
is thus given by VqJVq, = c ad /c en . 
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